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We investigate the distribution of the number of photons emitted by a single molecule undergoing 
a spectral diffusion process and interacting with a continuous wave laser field. The spectral diffusion 
is modeled based on a stochastic approach, in the spirit of the Anderson-Kubo line shape theory. 
Using a generating function formalism we solve the generalized optical Bloch equations, and obtain 
an exact analytical formula for the line shape and Mandel's Q parameter. The line shape exhibits well 
known behaviors, including motional narrowing when the stochastic modulation is fast, and power 
broadening. The Mandel parameter, describing the line shape fluctuations, exhibits a transition 
from a Quantum sub-Poissonian behavior in the fast modulation limit, to a classical super- Poissonian 
behavior found in the slow modulation limit. Our result is applicable for weak and strong laser field, 
namely for arbitrary Rabi frequency. We show how to choose the Rabi frequency in such a way that 
the Quantum sub-Poissonian nature of the emission process becomes strongest. A lower bound on 
Q is found, and simple limiting behaviors are investigated. A non-trivial behavior is obtained in the 
intermediate modulation limit, when the time scales for spectral diffusion and the life time of the 
excited state, become similar. A comparison is made between our results, and previous ones derived 
based on the semi-classical generalized Wiener-Khintchine formula. 
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I. INTRODUCTION 



Physical, Chemical, and Biological systems are inves- 
tigated in many laboratories using single molecule spec- 
troscopy The investigation of the distribution of 
the number of photons emitted from a single molecule 
source is the topic of extensive theoretical research e.g. 

S!SS!ffS!HmmE0andrj3fora 

review. Since optical properties of single molecules are 
usually very sensitive to dynamics and statics of their 
environment, and since the technique removes the many 
particle averaging found in conventional measurement 
techniques, single molecule spectroscopy reveals interest- 
ing fluctuation phenomena. An important mechanism 
responsible for the fluctuations in the number of photons 
emitted from a single molecule source is spectral diffusion 
e-g- 0, El E, EH E3- 111 many cases the absorption 
frequency of the molecule will randomly change due to 
different types of interactions between the molecule and 
its environment (e.g. EH HJ IH |H H |H H and 
Ref. therein). For example for single molecules embed- 
ded in low temperature glasses, flipping two level systems 
embedded in the glassy environment, induce stochastic 
spectral jumps in the absorption frequency of the single 
molecule under investigation 0, H3- In this way 
the molecule may come in and out of resonance with the 
continuous wave laser field with which it is interacting. 

Obviously a second mechanism responsible for fluctu- 
ations of photon counts is the quantum behavior of the 
spontaneous emission process [23, [3(J • In his fundamen- 
tal work Mandel [3l| showed that a single atom in the 
process of resonance fluorescence, in the absence of spec- 
tral diffusion, exhibits sub-Poissonian photon statistics 
[32j| . Photon statistics is characterized by Mandel's Q 



parameter 
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where N is the number of emitted photons within a 
certain time interval. The case Q < is called sub- 
Poissonian behavior, while Q > is called super- 
Poissonian. Sub-Poissonian statistics has no classical 
analog 31] . Briefly, the effect is related to anti-bunching 
of photons emitted from a single source and to Rabi- 
oscillations of the excited state population which fa- 
vors an emission process with some periodicity in time 
(see details below). Sub-Poissonian photon statistics 
and photon anti-bunching were measured in several 
sing le molecule, and sing le quantum dots experiments 
[lillllliSElllilla. While sub-Poissonian statis- 
tics is well understood in the context of resonance fluo- 
rescence of an isolated electronic transition of a simple 
atom in the gas phase, our theoretical understanding of 
sub-Poissonian statistics for a molecule embedded in a 
fluctuating condensed phase environment is still in its 
infant stages. 

In this paper we obtain an exact analytical expression 
for the Q parameter in the long time limit, for a sin- 
gle molecule undergoing a stochastic spectral diffusion 
process. To obtain the exact solution we use the Zheng- 
Brown generatin g fu nction method for single molecule 
photon statistics |40l liH l42j. For the spectral diffusion 
we use a simple stochastic approach, in the spirit of the 
Kubo-Andersen line shape theory [43l 0] . The model 
exhibits generic behaviors of line shapes of molecules 
embedded in a condensed phase environment, e.g. mo- 
tional narrowing when the stochastic fluctuations are 
fast, power broadening etc. We show that the Q pa- 
rameter exhibits rich types of behaviors, in particular 
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it reveals the quantum nature of the emission process 
in the sub-Poissonian regime, while the corresponding 
model line-shape exhibits a classical behavior. A brief 
summary of our results was recently published |45j . 

Our analytical expressions for Q classify the transi- 
tions between sub and super Poissonian statistics. They 
give the conditions on the spectral diffusion time scale 
for sub-Poissonian behavior. Motional narrowing type of 
effect is revealed also for the Q parameter. Our exact 
result is valid for weak and strong excitation (i.e arbi- 
trary Rabi frequency). It yields the lower bound on Q. 
The solution shows how we may choose the Rabi fre- 
quency so that the quantum nature of the photon emis- 
sion process becomes larger, namely how to minimize Q 
in the Sub-Poissonian regime. This is important for the 
efficient detection of quantum effects in single molecule 
spectroscopy, since choosing too small or too large val- 
ues of the Rabi frequency results in very small and hence 
undetectable values of Q. 

Finally our exact result is used to test the domain of va- 
lidity of the generalized Wiener Khintchine which yields 
Q in terms of a Fourier transform of a three time dipole 
correlation (as well known the Wiener Khintchine theo- 
rem yields the line shape in terms of a one time dipole 
correlation function). The theorem [2^, [2|| is based on 
the semi-classical theory of interaction of light with mat- 
ter, and on linear response theory (i.e., weak Rabi fre- 
quency), it yields Q > 0. As pointed out in [IE EE EH 
such a behavior is expected to be valid only for slow 
enough spectral diffusion processes. 



II. INTRODUCTION TO SUB-POISSONIAN 
STATISTICS 

We briefly explain some of the main ideas of sub- 
Poissonian statistics. The general idea is that the pho- 
tons emitted from a single particle, e.g. a molecule, a 
nano-crystal or atom are correlated in time. Consider 
first a hypothetical molecule, interacting with an ex- 
citing laser field, which emits photons with a constant 
time interval r between successive emission events. Then 
N = t/r, N 2 = N , and hence Q = — 1. Due to quan- 
tum uncertainty the photon emission process is always 
random and therefore — 1 < Q. Sub-Poissonian behav- 
ior where — 1 < Q < implies that the stream of pho- 
tons emitted from a single source maintain correlations 
in their arrival times to a detector. 

Usually when many molecules interact with a contin- 
uous wave laser the emission events are not correlated, 
and the fluorescence exhibits Poissonian statistics Q = 0. 
In contrast a single molecule, once it emits a photon, is 
collapsed to its ground state. Hence immediately after an 
emission event the molecule cannot emit a second time 
(it has to be re-excited by the laser). Hence successive 
photons emitted from a single molecule, seem to repel 
each other on the time axis, a non-Poissonian behavior. 
This well known effect is called anti-bunching 0, 



which is related to sub-Poissonian statistics. 

A second effect related to sub-Poissonian behavior are 
Rabi-oscillations. Consider a simple atom in the process 
of resonance fluorescence. When the electronic transition 
(frequency ojq) is in resonance with a continuous wave 
laser field (frequency ojl) the electronic transition can be 
approximated by a two level system. First let us men- 
tally switch off the spontaneous emission, i.e. set the in- 
verse life time of the transition T = 0. For zero detuning 
lol = ujq the transition will exhibit well known Rabi os- 
cillations: the population of the excited state will behave 
like p ee = sin 2 (Ot/2). Since the population in the ex- 
cited state attains its maximum (minimum) periodically, 
also the emission times of successive photons maintain 
certain degree of periodicity in time, which implies sub- 
Poissonian statistics. Mandel showed |3l| , that for a two 
level atom in the process of resonance fluorescence 

Q = 2- 2 

(r 2 + 2fF) 2 

When £1 <C T we have Q — > since then successive pho- 
ton emission times are not correlated, because the time 
between successive emissions becomes very large. While 
when !l>r the excited state is populated swiftly, and 
only the finite spontaneous emission rate delays the emis- 
sion, hence Q — > also in this case. Using Eq. the 
lower bound Q > —3/4 is easily obtained, and the mini- 
mum Q mm = —3/4 is obtained when ^ mm = T/y/2. 

III. MODEL AND GENERATING FUNCTION 
FORMALISM 

Let N be the random number of photons emitted by a 
single molecule source in a time interval (0, t), and Pjv(i) 
is the probability of N emission events. The information 
about the photon statistics is contained in the moment 
generating function |40| 

oo 

2y(s) = J2 sNp N(t) (3) 

which yields the moments of N 

N(f) = 2y'(i) jv5(i) = 2y" (i) + 2y (i), (4) 

with which the Q parameter can in principle be obtained. 
In Eq. (0J, and in what follows, we use the notation 

^-g( S )\ s=1 ^g'(l) ( 5 ) 

and similarly for second order derivatives with respect 
to s. The over-line in Eq. (0J describes an average over 
the process of photon emission, later we will consider 
a second type of average with respect to the spectral 
diffusion process, which we will denote with (...). 

The equations of motion for the generating function 
was given in 0] an d are called generalized optical Bloch 
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equation. For a chromophore with single excited and 
ground states, and interacting with a continuous wave 
laser field 

U (s) = -ZU (s) + 5 L (t)V (s) 

V (s) = -5 L {t)U (s) - EV (a) - ftW (s) 

W (s) = fiV (s) - § (1 + s) W (s) - | (1 + s) ^ (s) 

3>( S ) = -E(l- S )W( S )-E(l- S )y(s). 

(6) 

These equations are exact within the rotating wave ap- 
proximation and optical Bloch equation formalism. They 
yield the same type of information on photon statistics 
contained in the quantum jump approach to quantum op- 
tics which is used in quantum Monte Carlo simulations 
[30l l48| . In Eq. © T is the spontaneous emission rate 
of the electronic transition and Q is the Rabi frequency. 
The time evolving detuning is 



SL(t) — ujl - luq - Aw(t), 



(J) 



where ujl (ujq) is the laser frequency (the molecule's bare 
frequency), and Au(t) is the stochastic spectral diffusion 
process. In Eq. JBJ it is assumed that the molecule in 
its excited and ground state have no permanent dipole 
moments, hence the system is described only by the tran- 
sition dipole moment via the Rabi frequency. 

The physical meaning of U(s), V(s), and W(s) and 
their relation to the standard Bloch equation was given 
in ^fji some discussion on this issue will follow Eq. (|10|l . 
For related work on the foundations of these equations 
see 0, and Ref. therein. Note that when s — > 1 
the damping terms in Eq. © become small [i.e. the 
(1 — s)r/2 terms], hence relaxation of the generalized 
Bloch equations in the important limit of s — > 1 is slow. 

In what follows we will consider the moments N(t), 
N 2 (t). For this aim it is useful to derive equations of 
motion for the vector z = 

Mi) 1 v(i) 1 w(i)j(i) 1 w{i) 1 v(i) ) W(i) ) y(i) ) /(i)}. 

(8) 

Taking the first and the second derivative of Eq. © with 
respect to s and setting s = 1, we find 



i = M(t)z 
where M(t) is a 9 x 9 matrix M(t) = 



(9) 
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(10) 

The first three lines of M(t) describe the evolution of 
U{\), V(l), W(l), these are the standard optical Bloch 



equations in the rotating wave approximation [23 • These 
equations yield W(l) which in turn gives the mean num- 
ber of photons using Eq. Q 



N(t) = r 



(ii) 



and us- usual W(l) + \ is the population in the excited 

state. The fourth line of M(t) is zero, it yields 3^(1) = 0, 
this equation describes the normalization condition of the 
problem namely 3^(1) = 1/2 for all times t [to see this 
use Eq. © and J2n=o ^ n M = 1] • ^ ne evolution of 
the other terms W(l), V'(l), W'(l), y" (1) are of current 
interest since they describe the fluctuation of the photon 
emission process. In particular using Eqs. (|4I1U|) 



dt 



N(t) [N(t)-1] 



N (t) + 2W'(1) 



(12) 



Solutions of time dependent equations like Eq. © 
are generally difficult to obtain, a formal solution is 
given in terms of the time ordering operator T, z{t) = 
Texp[/*M(i)dt]z(0). 

Eq. jSj yields a general method for the calculation of 
Q for a single molecule undergoing a spectral diffusion 
process. The aim of this paper is to obtain an exact 
solution of the problem for an important stochastic pro- 
cess used by Kubo and Anderson [43L I53L l54| to investi- 
gate characteristic behaviors of line shapes. We assume 
Aco(t) — vh{t) where v describe frequency shifts, and 
h(t) describes a random telegraph process: h(t) = 1 or 
h(t) = —1. The transition rate between state up (+) and 
state down (-) and vice versa is R. This dichotomic pro- 
cess is sometimes called the Kubo- Anderson process [Sij . 
It was used to describe generic behaviors of line shapes 
here our aim is to calculate Q describing 
the line shape fluctuations. 

We use Burshtein's method [fill fs3 | of marginal aver- 
ages, to solve the stochastic differential matrix equation 
©. The method yields the average (z) with respect to 
the stochastic process. We will calculate (z) in the limit 
of long times, and then obtain the steady state behavior 
of the line shape 



(13) 



and the Q parameter. Let (z)± be the average of z(t) 
under the condition that at time t the value of h(t) = 
±1 respectively. (z)± are called marginal averages, the 
complete average is (z) = (z) + + (z)— . The equation of 
motion for the marginal averages is an 18 x 18 matrix 
equation 



(z) 



M + - RI RI 
RI M_ - RI 



(14) 



In Eq. i|14|) the matrix M± is identical to matrix M in 
Eq. l(TU|) when <5l(£) is replaced by 6^ = ujl — ujq T v, 
and J is a 9 x 9 identity matrix. In the next subsection 
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we obtain the long time solution of Eq. (|14f> . the reader 
not interested in the mathematical details may skip to 
subsection IIV Al where the solution for the line shape 
and Q is presented. 



IV. MATHEMATICAL DERIVATION OF 
EXACT SOLUTION 

In three main steps, we now find the long time behav- 
ior of the marginal averages, with which the line and Q 
are then obtained. 

1. As mentioned, from normalization condition we have 
y(l) = 1/2 for all times. Eq. (|14|l yields the marginal 
averages (3^(1)) + = (3^(1))- = 1/4, in the steady state. 
Inserting these identities in Eq. (|14fl we obtain an equa- 
tion of motion for the vector 

j, = ({W(l) + ,V(l)+,>V(l) + ,W(l)_,V(l)_,W(l)_}) 

(15) 



where A 
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(17) 



where A 3 § = ^4 63 1 , and A^ 1 is the ij matrix element of 
A- 1 , ij = 1,- • - ,6. We note that (W ss (l))± yields the 
steady state marginal averages of the population differ- 
ence between the excited and ground state. From Eq. 
|@J| we see that we need 2(3/(1)) to obtain the average 
number of photon emissions (N). We use Eq. (|14fl and 
show 

(y ss (i)) = <y ss (i))+ + (y ss (i))- = 



l^W ss (l)) + + {W ss (l))- + \)l (20! 

and 



<y-(i» + -<y-(i»- - r « w "")> : - <»"•"»-» . 

The average number of emitted photon, in the long time 
limit is lUEnj 



(N) = r ( (W ss (1))+ + (W ss (1))_ + - ) t, (22) 



namely (TV) is proportional to the steady state occupa- 
tion in the excited state. The line shape is 



I+ + I- 



and b = (0, 0, -r/4, 0, 0, -T/4). In the long time limit 
the solution of Eq. flfi)) reaches a steady state (ss) given 
by 

y" = -A-%, (18) 
and A -1 is the inverse of A. From Eq. I |18| l we find 



(W"{1))- = ~(A£ + A£), (19) 



where 



I± = 2T 



(w ss (i))± + i 



(23) 



(24) 



2. We use the solutions obtained in previous 
step to obtain inhomogeneous equations for x = 
({U'(1) + ,V'(1) + ,W'(1) + ,U'(1)-,V'(1UW'(1)-}), 

x = Ax + b(t) (25) 

where and b(t) = (0, 0, b+(t), 0, 0, &_(*)) with 



b±(t) = -I - ^<W SS (1)) ± - ^ {1 + 2 l(W ss (l)) + + (W ss (l))_]} ± ^ [<W SS (1))_ - (W ss (l)) + ] . (26) 
In the long time limit we obtain 

x(t)~A- 1 c + (A- 1 t + A- 1 A- 1 )c l (27) 

where cq and c\ are column vectors 



co = (0, 0, c+, 0, 0, c_) , Cl = (o, 0, 0, 0, ^p^) 



(28) 
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with 



C±= 2 



( W »(l)) ± + _ 



T^[(W SS (1))--(W-(1)) + ] 



we therefore obtain 



(W " (1)> + = A^c + + A^c- + I( WL )(W" (l))+t + [(^ 1 ^ 1 ) 33 + (^^J 



(29) 



(30) 



(W ss (1))_ = A£c+ + A£c- + I(u L )(W ss (l))_t + 



r/ (u £ ) 



[(^- 1 )« + (^-^- l ) M ] 



(31) 



3. From Eq. Q, (N 2 (t)} = (2y"(l)) + (2y'(l)). In steady state we have 

(y" ss (i)) = (y / "(i))+ + (y , "(i)>_. 

From Eq. (|14|l one can show that in the long time limit 

(y" ss (i)> ~ 4 [(w ss (i)> +c+ + (w ss (i))_ c _] t+ 



(32) 



r 2 /k) 



[(A-'A-% 3 + (A-'A-\ 6 + (A~iA-i) 63 + (A-'A-^J t 



r/K) 



Finally we obtain the Q parameter using: 



(W ss (l))+ + (W SS (1))- 



(y"(i) ss ) -2( y(i) ss ) 2 
(y(i) s 



(33) 



(34) 



Using Eqs. (|20I33I34|) we obtain the main result of this manuscript 

">2 T^2 y 



i— 3,6 j— 3,6 



ii/(wL) 



^{(W-(l)) fc [l + 4(W ss (l)) fc ]}, (35) 



fe=± 



which is valid when measurement time t — * oo. The Q parameter in Eq. I|35(l is expressed in terms of A -1 . To obtain 
the solution in terms of the original parameters of the problem i?, is, lul, luo, ft, T, we found analytical expressions for 
A^ 1 using Mathematica. The formula for the Q parameter i s given in the following subsection. 



A. Exact Solution 

Without loss of generality we set ojq = 0, hence u>l is 
the detuning. We find 



Numerator [C 
Denominator [ 



Numerator [Q] = 



(36) 



-2Tn 2 (R(T + 4i?)(3r 3 - 4Tv 2 + 24T 2 R + \Qv 2 R - 

48ri? 2 )(r 3 + er 2 i? + 8v 2 r + 2V{2v 2 + ft 2 + m 2 )) 3 

8(4r n i? - AY w {v 2 - 21R 2 ) - 3072r// 4 i? 5 (4i/ 2 - ft 2 
16i? 2 ) - 4096z/ 4 i? 6 (2t/ 2 - ft 2 - 4i? 2 ) + 15T 9 R(Q 2 
48R 2 ) - 8T 8 (4v 4 - 27ft 2 i? 2 - 392R 4 + 2^ 2 (ft 2 - 48i? 2 )) - 



4r 5 i?(128i/ 6 + ft 6 + 48ft 4 i? 2 + 1344ft 2 i? 4 + 6400i? 6 - 
12z/ 4 (llft 2 - 272E 2 ) + 4z; 2 (9ft 4 - 302ft 2 i? 2 - 5072i? 4 )) - 
64r 3 i? 3 (48i/ 6 + 3v 4 (llfl 2 - 448i? 2 ) + 3(ft 2 + 4i? 2 ) 2 (ft 2 - 
16i? 2 ) + 2^ 2 (5ft 4 - 80ft 2 i? 2 - 544E 4 )) - 256r 2 i? 4 (28^ 6 - 
*/ 4 (3ft 2 - 332i? 2 ) + (ft 2 + 4i? 2 ) 3 + 2i/ 2 (ft 4 - 4ft 2 i? 2 - 
32i? 4 )) + 4T 7 R(-48v 4 - 2i/ 2 (9ft 2 - 952i? 2 ) + 3(ft 4 - 
89ft 2 i? 2 + 544i? 4 ) ) - 8r 6 (8j/ 6 + 8u 4 (ft 2 - 17i? 2 ) + 2v 2 (ft 4 - 
38ft 2 i? 2 - 2116i? 4 ) - 3(3ft 4 i? 2 + 58ft 2 i? 4 + 40i? 6 )) - 
16r 4 i? 2 (96^ 6 + 96^ 4 (ft 2 - 31i? 2 ) + i/ 2 (26ft 4 - 720ft 2 i? 2 - 
6656i? 4 ) + 3(ft 6 + 44ft 4 i? 2 + 448ft 2 i? 4 + 1152i? 6 )))w| - 
32(r+2i?)(3r 8 i?+512r^ 4 i? 4 +512j/ 4 J R 5 -8r 7 (i/ 2 -3i? 2 )- 
r 6 (-32z/ 2 i? + 3i?(ft 2 -4i? 2 ))-32r 2 i? 3 (-14z/ 4 + 6^ 2 ft 2 - 
3(ft 2 + 4i? 2 ) 2 ) + 2r 5 (16^ 4 - 8^ 2 (ft 2 - 3i? 2 ) - 3i? 2 (5ft 2 - 
104i? 2 )) - 8r 3 i? 2 (-56^ 4 + v 2 (29ft 2 - 40i? 2 ) + 6(ft 4 - 
20ft 2 i? 2 + 64i? 4 )) - 2r 4 i?(-104i/ 4 + lfV 2 (5ft 2 - 16i? 2 ) - 
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3(fi 4 + mn 2 R 2 + 368i? 4 )))w| - i28r 2 (r + 2R) 2 {at 2 (v 2 - 

3R 2 ) + TR{l&v 2 + 3ft 2 + A8R 2 ) + AR 2 {Av 2 + 3(Q, 2 - 
AR 2 )))lu 6 l - 256r 2 i?(r + 2R) a Lu s L ) 

Denominator [Q] = 

R((T + AR)(T 3 + 6r 2 i? + 8v 2 R + 2T{2v 2 + il 2 + AR 2 )) - 

4r(r + 2R)lo 2 l ) (r 6 + ior 5 i? + §Av 2 si 2 r 2 + ay 4 (2v 2 + n 2 ■ 

8R 2 ) + AT 3 R(12v 2 + 7fl 2 + 8i? 2 ) + \&TR{2v 4 + 3v 2 n 2 - 

n 4 + in 2 R 2 ) + 4r 2 {Au A + n 4 + mn 2 R 2 + a v 2 {n 2 +m 2 ))- 
8r(r 3 + 6r 2 i? - 8^ 2 i? + 6Q 2 r + 16R 3 + 2T(-2v 2 + n 2 - 
8i? 2 ))^ 2 + i6r(r + 2i?K) 2 . 

The line shape is 

Numerator[7 (ojl)) 



,x 10 



Denominator [I (u>l)] 



(37) 



Numerator [I (cul)] 



rn 2 ((r + 4i?)(r 3 + 6r 2 i?+ 8^ 2 i?+ 2r(2^ 2 + n 2 + ar 2 )) 
4r(r + 2i?K), 

Denominator [I (u>l)] = 

r 6 + ior 5 i? + 64^ 2 ^ 2 i? 2 + AT i (2v 2 + n 2 + 8R 2 ) 

AT 3 R(12v 2 + 7Q 2 + 8R 2 ) + lQTR(2v 4 + 3v 2 Q 2 + fl 4 - 
4ft 2 i? 2 ) + AT 2 (Av 4 + fl 4 + l6fl 2 R 2 + Au 2 {VL 2 + AR 2 )) - 

8r(r 3 + 6r 2 i? - 8^ 2 i? + 6Q 2 r + 16R 3 + 2T{-2v 2 + n 2 

8R 2 ))uj 2 l + 16T(r + 2R)u>i. 




0> L - G> Q ( MHz ) 



1. Exact Solution for Zero Detuning 



When the detuning is zero we find 



2m 2 



Q(w L =0) = — 



4^ 2 (-r + AR) + 3T (r + ary 



[ATv 2 + (r 2 + 29?) (r + AR)Y 

(38) 

Using Eq. (J3HJ the lower bound Q (uj l = 0) > -3/4 is 
obtained. The absolute minimum of Q, i.e. Q = —3/4 is 
found when v = 0, Q — r/%/2, or when R — > 00, Q = 
Yj\[2. Namely the absolute minimum is found for a 
molecule whose absorption frequency is fixed, or for a 
very fast spectral modulation. 

Remark: Eq. I|38|) indicates a transition from sub- 
Poissonian statistics (Q < 0) to super-Poissonian statis- 
tics (Q > 0) when 4i/ 2 = 3(r + AR) 2 /(T - AR). If 
R > r/4, i.e. if the bath is fast compared with the 
radiative life time, we find sub-Poissonian behavior for 
all values of v and f2. 



V. THE PHYSICAL BEHAVIORS OF THE 
EXACT SOLUTION IN LIMITING CASES 

The Q parameter is a function of two control param- 
eters ojl and O, and three model parameters r, R, v. 
The classification of different types of physical behaviors, 
based on the relative magnitude of the parameters is in- 
vestigated in this section. The limiting behaviors of Q 
are obtained from the exact solution using Mathematica. 



FIG. 1: The line shape I(lol — <^>o) (a) and the Q(cjl — Wo) 
parameter (b) are calculated with the exact solution Eqs. I|M7[I 
and Eq. I36H respectively. In this slow modulation strong 
coupling limit R « T ~ Q « v, both Q and I exhibit 
splitting behavior with two peaks centered at ujl — uio = ±f. 
The Q parameter exhibits super-Poissonian statistics Q > 0. 
The * asterisks are approximate results Eq. I|39|l for the line 
shape, and Eq. H41H for Q. Parameters are T — 40 MHz, 
v = 5r, n = T/V2, and R = 1 Hz (solid curve), 5 Hz (dashed- 
dot curve), and 25 Hz (dashed curve). 



A. Slow Modulation Regime: R « v, T, H 

In the slow modulation regime, the bath fluctuation 
process (R) is very slow compared with the radiative de- 
cay rate (r), frequency fluctuation amplitude (v) and 
the Rabi frequency (Q). This case is similar to situations 
in many single molecule experiments, for example single 
molecules in low temperature glasses. 

The exact solution can be simplified in the limit R — > 0, 
using Eq. JjSJ 



lim I (lu l ) 



1+ (uj l ) + I- {uj l ) 



and 



ra 2 

(r 2 + 2r» 2 + 4 (uj l t- vf) 



(39) 



(40) 



7 



The line is a sum of two Lorentzians centered on ±i/, 
namely it exhibits splitting behavior when v » r,ft. 
Using Eq. (|3"r))> we find a simple super-Poissonian behav- 
ior for Q in the slow modulation limit R — > 

Q-Qslow^ (/+ 4 flj~ )2 > ( 41 ) 

and since R — > 0, Q may become very large (e.g. Q — 
5 * 10 6 in Fig. [TJ. 

A simple picture can be used to understand these re- 
sults. In the slow modulation limit the molecule jumps 
between two states + and — , the time between succes- 
sive jumps is very long, in such a way that many photons 
are emitted between jump events. In each of these two 
states the molecule emits photons at a rate I±{u>l)i Eq. 
l|4Up. These rates are determined by the familiar steady 
state solutions of the optical Bloch equation, for a two 
level atom with the absorption frequency luq ± v fixed in 
time |28j . In this slow limit the random number of emit- 
ted photons, in time interval (0,i), is iV g j ow = L I(t)dt, 
and I(t) is a stochastic intensity that jumps between two 
states I±(u>l) with the rate R. Using this simple random 
walk picture it is straightforward to derive Eqs. I|39l I41|l . 
For mathematical details see Ref. [2^ who considered 
a similar slow modulation limit which is valid only for 
weak Rabi frequency. 




-200 -150 -100 -50 50 100 150 200 
to L - co Q ( MHz ) 



i. r « r ~ << v 

Within the slow modulation limit we distinguish 
between two cases. The case R << T ~ ft << v is 
called the slow modulation strong coupling limit. In 
this case the line and Q have two well separated peaks 
and the broadening of the two peaks due to the finite 
life time, and power of the laser field is small compared 
with the frequency shifts (see Fig. QJ. From Fig. ^ and 
as expected from Eqs. l|3lfl l4~T)l . Q decreases when R is 
increased, while I(u)l) is independent of R. Thus it is 
Q not I{<j->l) that yields information on the dynamics. 
In Fig. ^the agreement between the exact solution and 
the approximation Eq. lj4"T|) is good. 



2. r« v <<r ~n 

The limit i?<<^<<r~ftis called the weak cou- 
pling slow modulation limit. In this case the two peaks of 
the line, discussed in the previous sub-section, are over- 
lapping and the line is approximated by 



Tft 2 



r 2 + 2ft 2 



4< 



(42) 



This result is exact when v = 0, for arbitrary R, 

The behavior of Q is demonstrated in Fig. [3 where 
we observe both super-Poissonian and sub-Poissonian be- 
haviors. In the slow modulation weak coupling limit, we 



FIG. 2: Same as Fig. 0for the slow modulation weak coupling 
case R « v « F ~ Q. We see that the line / is Lorentzian 
in shape while the Q parameter exhibits splitting. And unlike 
the line shape Q depends on R in this limit. In the vicinity of 
zero detuning Q exhibits a sub-Poissonian behavior (Q < 0). 
Far from zero detuning the slow spectral diffusion process 
controls the behavior of Q and then Q > 0. The * asterisk 
are the approximate Eqs. for Q Eq. 145H and I Eq. 1421 . 
The parameters are F = 40 MHz, v = T/10, SI = r/V2, and 
R = 172500 (solid curve), T/500 (dashed-dot curve), and 
I 1 / 100 (dashed curve). 



must distinguish between the cases of large and small de- 
tuning. First note that according to Eq. I|41(l when the 
detuning is zero we find Q — 0, namely the leading 
term in our asymptotic expansion vanishes. We must 
therefore consider the higher order terms in our asymp- 
totic expansion of Eq. (|36|l and we find 



Q ~ Qslow + Qm 



where 



2ft 2 (3r 2 - 4w 2 ) 

(r 2 + 2ft 2 + 4w 2 ) 2 ' 



(43) 



(44) 



Eq. I|43|) has a simple meaning, the first term is a con- 
tribution to Q from spectral diffusion, which is identical 
to Eq. (J2J. The second term Qy^ is identical to the 
result obtained by Mandel, for the Q parameter in the 
absence of spectral diffusion ; 3I], and Qm < provided 
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that the detuning is not too large. The second term is 
dominating over the first when the detuning is small, and 
for zero detuning we obtain in Eq. H43|l sub-Poissonian 
statistics. More explicitly, we Taylor expand (|40|l using v 
as a small parameter, and obtain for the slow modulation 
weak coupling limit 



Q 



MTti 1 uj 2 T v 2 



2ft 2 



( 3 r 2 _ 4^1) 



r (r 2 + 2ft 2 



(r 2 + 2ft 2 



(45) 



One may say that for zero detuning, the molecule behaves 
as if its absorption frequency is fixed. 

To conclude, we see that in the slow modulation limit 
Q is a sum of two additive contributions: i) a part related 
to spectral diffusion Q s i ow and ii) and a part related to 
quantum fluctuations , i.e . Qm- Such a behavior was very 
recently discussed in ( see also [111 f° r related discus- 
sion). The quantum fluctuations are however bounded 
from above and below —3/4 < Qm < 0, while the con- 
tribution from spectral diffusion is not < Q s \ ow < oo. 
Hence detection of the quantum fluctuations is possible 
only when Q s \ ow is small, which for our case implies zero 
detuning and weak coupling limit. 



B. Fast Modulation Regime: R » v, Y, Q 




-400 -300 -200 -100 100 200 300 400 
co L - (o Q ( MHz ) 



We now consider the fast modulation limit. If we take 
R — > oo keeping ft, T, and v fixed we find from Eq. 136|) 



Jim Q = Qm 



(46) 



given in Eq. I|44|). Hence in this limit we find a sub- 
Poissonian behavior, provided that the detuning ljl is 
not too large. The line shape is identical to the expres- 
sion on the right hand side of Eq. (|42(l . This behavior 
is expected, when the spectral diffusion is very fast the 
emitting single molecule cannot respond to the stochastic 
fluctuations. 

A more interesting case is to let R — > oo and v — > oo 
but keep, 



rSD= R 



(47) 



finite. We call this limit the fast modulation limit, using 
Eq. (|3*7I) the line shape is 



J fast (^l) = 



(r + r SD )ft 2 



(r SD + iy + 2 (i + r SD /r) ft 2 + 4w 2 



and when the Rabi frequency is small 

(r 



1 fast 



r SD ) ft 2 



(r SD + ry + M 



(48) 



(49) 



In this limit we have the well known effect of motional 
narrowing: the width of the line is determine by T + 



FIG. 3: Same as Fig. 0for fast modulation R » v >> 
Q, — r/\/2. We see quantum fluctuations namely Q < 0. 
Both for the line shape and for Q we have motional narrowing 
behavior, as the spectral diffusion process becomes faster Q 
and / narrow. When R — > oo, Q approaches Mandel's formula 
Eq. H44fl . I n the same limit the width of the line shape is 
determined only by its natural width, namely by the inverse 
radiative life time of the excited state (provided that the Rabi 
frequency is weak). The * asterisk are the approximations for 
Q Eq. dSm and I{ujl) Eq. The parameters are r = 40 

MHz, v = 5r, Q = r/VI, and R = 5T (dashed curve), 25r 
(dashed-dot curve), and 125r (solid curve). 



Tsd and as the process becomes faster the line becomes 
narrower, namely Tsd decreases when R is increased. 
The Q parameter is obtained from Eq. H3(jfl 



ast 



2rft 2 [3r 3 



5rr 2 D 



SD 



7r 2 



[r 3 + rr SD 2 + 2rft 2 + 2r SD (r 2 + ft 2 ) 



(50) 

Thus in the fast modulation limit the photon statistics 
is sub-Poissonian provided that the detuning is not too 
large. When Tsd — ► the result for Q reduces to Man- 
del's result Eq. (|4*4f> . In Fig. |3| we show the line shape 
and the Q parameter, for three values of the jump rate 
(R) in the fast modulation regime while v, ft, T are kept 
fix. We see that as the stochastic spectral diffusion pro- 
cess gets faster, both the line shape and the Q parameter 
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FIG. 4: We show the transition of the Q parameter from sub- 
Poissonian behavior, to super Poissonian behavior as the rate 
R is changed. Note that Q exhibits a non-trivial behavior: it 
may have four peaks. Such a behavior is limited to a small 
regime of parameters when R ~ T ~ Q and was not observed 
in other limits of the problem. The parameters are P = 40 
MHz, v = 5r, Q = V/V2, and R = T/6 (solid curve), T/4 
(dashed-dot curve), and P/2 (dashed curve). 



J 



become narrow. Thus both I(u>l) and Q exhibit a mo- 
tional narrowing effect. Also, as the stochastic process 
gets faster, a stronger quantum behavior is obtained, in 
the sense that the minimum of Q decreases. 



C. Strong coupling limit v >> R, P, 



To investigate the strong coupling limit wc consider 
the value of Q for u>l — and v >> T, R, O. From Eq. 
(|36J) we obtain 



lim Q u 



(r + 2R) Tfl 2 [-T 3 + 16TR 2 + 8R 3 - 2Tfl 2 + 2R (2T 2 + ft 2 )] 



2R [r 3 + 4ri? 2 + 2Tfl 2 + 2R (2F 2 
I 



(51) 



This equation exhibits both sub-Poissonian and super- 
Poissonian behaviors. When the process is very slow, 
namely R — > 0, we obtain 



lim Q u 



Tfl 2 



2R(T 2 + 2n 2 ) 



(52) 



a super-Poissonian behavior. In the intermediate modu- 
lation limit when R = T, we obtain 



lim Q LUL=L , = - 



8m 2 r 2 



2 (9T 2 + 4fi 2 



(53) 



a sub-Poissonian behavior. 
Q is small 



lim Q u 



When R — > oo wc find that 



n 2 

2RY' 



(54) 



Eqs. H52I54|I are valid only in the limit of v — > oo. How- 
ever behaviors similar to the predictions of these equa- 
tions are found also for finite values of v. On lul = v wc 
have three typical behaviors: (i) Q{lol = v) > when 
the process is slow, see Fig. HJ. (ii) Q{lol = v) < when 
r ~ R, see Fig. 0]for R — T/2, and (Hi) when R — > oo 
we find Q (lo l = v) -»• 0, Fig. 
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FIG. 5: Sub - Poissonian behaviors of Q for the case of strong 
coupling v >> R and intermediate modulation R ~ T ~ 
limit. Q exhibits two types of quantum sub-Poissonian 
statistics: i) splitting behavior and ii) motional narrowing 
behavior where Q has a single minimum. The parameters 
are: F = 40 MHz, v = 5r, Q = T/v^, and i? = IT (dashed 
curve), 3r (dashed-dot curve), and 5r (solid curve). 



D. Intermediate Modulation Limit R ~ V 

When R ~ T ~ £1 << v we obtain interesting behav- 
iors for Q. In Fig. 0]the Q parameter shows a transition 
from sub-Poissonian to super Poissonian, photon statis- 
tics. In this regime of parameters, the shape of Q when 
plotted as a function of u>l is very sensitive to the value 
of the control parameters e.g. in Fig. Q]we change R only 
moderately still we see very different types of behaviors 
for Q. For certain values of parameters Q attains more 
than two peaks (see Fig. 0Jfor R = T/6, T/4). In contrast 
exhibits a simple splitting behavior with two peaks 
on ±i/, which is similar to the slow modulation case. 

Besides the transition from sub to super Poissonian 
behavior, a second type of transition is observed as R 
is increased. In our problem we have two types of sub- 
Poissonian behavior. We noticed already that when the 
stochastic modulation becomes very fast, Q has one min- 
imum on zero detuning (see Fig. while when R ~ V, 
Q has two minima on u>l ± v (see Fig. [3]and R = T/2). 
The transition between these two types of sub-Poissonian 




(0 L - tO Q { MHz ) 

FIG. 6: We investigate the dependence the line shape and 
Q on the Rabi frequency fl. Power broadening of the line 
is observed, while for the Q parameter we have a turn-over 
behavior as explained in the text. The three values of the 
Rabi frequency are O = T/A (dashed curve), F/V2 (dashed- 
dot curve), and ST/2 (solid curve). And the minimum of Q, 
on zero detuning, is obtained for the intermediate value of 
n = r/v/2. The fixed parameters are r = 40 MHz, v = 5r, 
and R = WOT. 

behaviors is shown in Fig. 



VI. EXTREMUM OF Q 

We now investigate the dependence of Q on the ex- 
citation field. In Fig. we consider an example line 
shape and Q parameter, where we fix the model parame- 
ters v, i?, r and vary the Rabi frequency. For the line we 
see well known power broadening: as the Rabi frequency 
is increased the line becomes wider, and as expected the 
photon emission rate I{ujl) increases monotonically when 
fl is increased. For the Q parameter we have a turn-over 
behavior, as we increase f2 the value of Q on zero detun- 
ing decreases then increases. 

Generally this type of turn-over is expected, since as 
discussed in Sec. [H] Q = when ft — > oo or £1 — > 0. Thus 
there exists an optimal Rabi frequency which yields an 
extremum of Q. Obviously it is important to obtain the 
values of which yield the extremum of Q, since then 
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FIG. 7: The photon emission rate 1(H) (a) the Q(0.) param- 
eter (b) are calculated with the exact solution of Eq. 13 7H and 
Eq. J^J, in the fast modulation limit. The parameters are 
T = 40 MHz, v = 5I\ lu l = 0, Q = -> 7r, and R = 10F. We 
observe the saturation of the line as Q becomes large, while 
Q(£l) has a minimum. The value of Q which minimizes Q(Q) 
is of interest since it yields the strongest quantum fluctua- 
tions. The symbol * is the point (Sl mm , Qmin) calculated 
based on the approximation Eqs. (15715811 . 



the fluctuations are the largest. The extremum can be 
either a minimum or a maximum, as we shall show now. 



o 




FIG. 8: We demonstrate how to choose the Rabi frequency 
in order to obtain strong sub-Poissonian behavior. The Q(S1) 
parameter is calculated with the exact solution of Eq. 1361 . 
The detuning ujl — <^o is zero, and T = 40 MHz, v — F/10 
and R = T/100 corresponding to the slow modulation weak 
coupling limit. The minimum, the * in the figure, Q mm — 



-0.71 is found for £!_ 



minimum (Q 



and 



rum' ^min) 



0.72r, which is close to the global 

:(-3/4,TA/2). 



Uv 2 (-Y + AR) + 3r(r + 4i?) : 

4 (T + AR) (r 2 +4v 2 + ATR) 



(56) 



Eqs. (|55I56I> yield £l mm and Q mm in terms of v and 
R. Due to motional narrowing effect, for fast processes 
satisfying v << R, v and R are not easily obtained from 
experiment, while the parameter Tsd is in principle easy 
to obtain from the measurement of the line width. Using 
Eq. (|50|l . we find in the fast modulation limit and for 
zero detuning 



a 



i + r SD /r 



(57) 



A. Largest Quantum Fluctuations 

We now consider the quantum regime Q < 0. In Fig. 
Ulb) we demonstrate the turn-over behavior of Q(il) for 
an example where the stochastic fluctuations are fast. 
In this fast modulation case Q < hence Q(£l) has a 
minimum. For the same parameters the photon emission 
rate I(£l) saturates as is increased, and the emission 
rate is never faster than T [see Fig. 0a)]. 

Let ^ mm be the Rabi frequency which minimizes Q 
in the sub-Poissonian case Q < 0, and Q m i n the corre- 
sponding value of Q. Using Eq. I|38() we find for zero 
detuning, and for R > T/4 



'r 3 + 4r 2 i? + 4r^ 2 
2 (r + Ir) 



(55) 



(r SD + 3r) 
4(r 



SD) 



(58) 



These simple equations relate between the width of the 
line given in Eq. (|49J) and Q m [ n and ^ mm . From Eqs. 
(|57I58|I we see that when T SD « T, ft min = T/y/2 
When T S d >> T we find fL 



and Q 



mm 



-3/4. 



"mm 



V r s D T/2 and Q min = -1/4. 

Remark: When the detuning is not zero, we find using 
Eq. (EOJ 



' (i + r SD AT + 4i. 2 /r 2 
2(i + r SD /r) 



(59) 



Similar turn-over behaviors of Q(£l) are found also in 
other non fast parameter regimes. In Fig. [S] we show 
Q(Q) versus 51 for the slow modulation weak coupling 
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FIG. 9: Same as Fig. El for the intermediate modulation 
limit F ~ R. Unlike Figs. [7E1 the detuning is on ujl = v. 
Parameters are chosen as T = 40 MHz, v = 5r, SI = — > 20r, 
and R — 2T. The figure illustrates that the turn-over behavior 
of Q(fi) is generic, and not limited to the fast modulation 
limits Fig. |7|and slow modulation/weak coupling limit Fig. 

m 



limit R < v < r and for zero detuning. Fig. [SJshows that 
Q exhibits a minimum as function of SI, this minimum is 
found in the vicinity of f2 mm = T/y/2. Such a behavior is 
understood based on Eq. i|43|) , the spectral diffusion con- 
tribution for Q is not important at zero detuning, while 
the contribution of Qm yields ^ mm = T/ ^/2. To demon- 
strate that the turnover behavior of Q(Q) is generic, we 
consider also the intermediate modulation limit in Fig. [5] 
Here we choose the detuning according to ujl = v, since 
the Q parameter on zero detuning is relatively small (see 
Fig. 0). 



B. Maximum of Super-Poissonian Fluctuations 

In contrast to the behaviors in the quantum regime 
Q < 0, in the slow modulation limit where Q > 0, Q(S1) 
obtains a maximum, whose location is easy to calculate 
with Eq. H41f) . Such a behavior is demonstrated in Fig. 
IHJI for a case when the spectral shift v is not very large. 
If v >> SI, r then in the slow modulation limit 



Qslow 



rsi 2 



2R (T 2 + 2Q 2 



3rSl 2 
32Rv 2 



(60) 



when the detuning is equal to lu^ = v. In Eq. I|6(J|) 
the second term on the right hand side is supposed to 
be a correction to the first term, namely Q s i ow > 0. 
Let Sl m ax be the value of SI which maximizes Q in the 
super-Poissonian regime, and Qmax the corresponding 
maximum. This maximum always exists since as men- 
tioned Q = when CI — > or Q — > oo. Then using Eq. 
(H3J) 



max 



T(4V3f-3r) 



(61) 



FIG. 10: The crossover behavior of the Q (SI) parameter 
in the slow modulation limit. Q exhibits a super-Poissonian 
behavior, and Q (Q,) attains a maximum when the Rabi fre- 
quency is changed. Thus also for this classical type of behav- 
ior an extremum of the fluctuations is found for a particular 
value of the control parameter f2. The parameters are r = 40 
MHz, v — 5r, R = 500Hz and ujl — Wo = v. Using the ap- 
proximate Eqs. I|61I62|I we obtain the values Qmax — 16686 
f2max — 2.3r, the symbol * in the figure, which is in good 
agreement with the exact solution (the solid curve). 



which is independent of R and 



Qi 



r (3r 2 + 16^ 2 - 8V3IV) 



64i?^ 2 



(62) 



Note that when the frequency shifts are very large v — > oo 
we find using Eq. (|61|l ilmax - ¥ oo. Hence the value of 
Jl m ax may become very large and then in experiment it 
is impossible to reach fimax (e.g. v — \GHz). If we 
impose the condition SI << v we have 



*3slow 



Tfl 2 



2R{T 2 + 2tt 2 



(63) 



for the laser detuning ujl = v. Hence Q s i ow monoton- 
ically increases and eventually saturates, similar to the 
behavior of the average emission rate. 



VII. SUMMARY AND DISCUSSION 

The Q parameter yields informations not contained in 
the line shape. The most obvious is the transition from 
super (i.e. classical) to sub (quantum) Poissonian behav- 
ior. Such a quantum signature of the photon emission 
process is not obtained from the line shape. In compar- 
ison with the Q parameter of a single atomic transition, 
the Q parameter investigated here exhibits rich behav- 
iors. These include splitting, both in the sub and in the 
super Poissonian regime, a transition from a fast to a 
slow modulation limit, and effects related to motional 
narrowing. The most non-trivial behavior is obtained 
in the intermediate modulation limit when r ~ R ~ SI 
where Q attains more then two peaks. 



13 



Since Q contains the new information on single 
molecule experiments, namely information beyond the 
line shape, it is important to emphasize that Q attains an 
cxtrcmum for a particular value of the Rabi frequency. In 
particular in the sub-Poissonian regime Q(fl) has a min- 
imum. Hence we optimize the Rabi frequency in such 
a way that \Q\ is increased, e.g. we obtain £^ mm - In 
other words there exist an "ideal" choice of the Rabi fre- 
quency in single molecule experiments. In the quantum 
sub-Poissonian regime this optimal Rabi frequency can- 
not be considered weak, neither strong, hence perturba- 
tive approaches to single molecule spectroscopy are not 
likely to yield it. This is in complete contrast to most 
theories of line shapes which are based on the assump- 
tion of weak external fields, e.g. the Wiener-Khintchinc 
theorem and linear response theory. Single molecule the- 
ories should be able to predict the turn-over behavior of 
Q(f2) based on different models, since such a behavior 
is not expected to be limited to the model under inves- 
tigation. Of-course the exact solution presented in this 
manuscript is very valuable in this direction, since it pre- 
dicts precisely the details of this transition for the Kubo- 
Anderson stochastic process. 

It would be interesting to investigate further how gen- 
eral are our results. From line shape theory, we know that 
in the fast modulation limit, line shapes have Lorcntzian 
shapes under very general conditions. From experiment 
we know that motional narrowing effect, and Lorcntzian 
behavior of lines is wide spread. Thus at-least in this 
limit certain general features of line shapes, which are 
not sensitive to model assumptions are found. Similarly, 
we expect, that in the fast modulation limit, some of our 
results are general. For example the motional narrowing 
behavior of Q and its approach to Mandel's behavior Qm 
is likely to be general. It would be interesting to check if 
the relation between <5 mm and ^ m j n Eq. (|57I58(I . and 
the width of the line given by T$d and T is valid for 
other models, both stochastic and Hamiltonian. These 
simple equations are important since they yield the op- 
timal Rabi frequency ^ mm in terms of the width of the 
line shape, which in turn is easily determined in usual 
line shape measurement. 

In Ref. (|2{|,|26j) a semi-classical framework for the 
mathematical calculation of Q for single molecule spec- 
troscopy was investigated. The approach yields the Q 
parameter in terms of a Fourier transform of a three 
time dipole correlation function. As mentioned in the 
introduction, the approach in (H^Hil) is based on two 
main approximations (i) external fields are weak r2 — ^ 0, 
i.e. linear response theory, and the (ii) semi-classical ap- 
proach to photon counting statistics. The second as- 
sumption implies that Q > 0, and as pointed out in 
|lfil l2("il| such an approach is expected to be valid for 
slow processes. The approach is useful since most single 
molecule experiments report on slow fluctuations. The 
results obtained in this manuscript reduce to those in 
( (2^1, [2f|) in the limit of f2 — * 0, and in the slow modu- 
lation limit, as they should. The quantum behavior of Q 



becomes important when R ~ T or for faster processes. It 
is left for future work to construct a general quantum lin- 
ear response theory, based on the Eqs. of motional which 
would yield both super and sub-Poissonian statistics. Fi- 
nally, also the investigation of the time dependence of Q 
is timely. 

This work was supported by National Science Foun- 
dation award CHE-0344930. EB thanks the Complexity 
Center in Jerusalem for financial support. 
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